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Abstract 

Unit-vector fields n on a convex polyliedron P subject to tangent 
boundary conditions provide a simple model of nematic liquid crystals 
in prototype bistable displays. The equilibrium and metastable config- 
urations correspond to minimisers and local minimisers of the Dirichlet 
energy, and may be regarded as S -valued harmonic maps on P. We con- 
sider unit-vector fields which are continuous away from the vertices of P. 
A lower bound for the infimum Dirichlet energy for a given homotopy class 
is obtained as a sum of minimal connections between fractional defects at 
the vertices of P. In certain cases, this lower bound can be improved by 
incorporating certain nonabelian homotopy invariants. For a rectangular 
prism, upper bounds for the infimum Dirichlet energy are obtained from 
locally conformal solutions of the Euler-Lagrange equations, with the ratio 
of the upper and lower bounds bounded independently of homotopy type. 
However, since the homotopy classes are not weakly closed, the infimum 
may not be realised; the existence and regularity properties of continuous 
local minimisers of given homotopy type are open questions. Numerical 
results suggest that some homotopy classes always contain smooth min- 
imisers, while others may or may not depending on the geometry of P. 
Numerical results modelling a bistable device suggest that the observed 
nematic configurations may be distinguished topologically. 

This article appears as a chapter in "Analysis and Stochastics of Growth 
Processes and Interface Models", P Morters et al. eds., Oxford University 
Press 2008, http://www. oup. com/uk/catalogue/?ci=9780199239252 . 



1 Introduction 

Liquid crystals are intermediate phases of matter exhibiting partial ordering in 
the orientation and/or positions of their constituent particles. The constituents 
of nematic liquid crystals have a distinguished axis, and in the nematic phase 
these axes tend to align. The direction and degree of alignment can exhibit 
a rich variety of singularities. Standard references on liquid crystals include 

H [Ml 1121]. 
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Figure 1: The PABN cell (from [5]) 



The nematic phase is optically birefringent (light propagation is polarisation- 
dependent) . This, together with the fact that nematic ordering can be modified 
by external electric and magnetic fields, has led to a wide range of display appli- 
cations. Most present-day liquid crystal displays (eg twisted nematic) are based 
on monostable cells, where, in the absence of external fields, the orientation 
assumes a single (spatially varying) equilibrium configuration which is effec- 
tively transparent to incident polarised light. To produce and maintain optical 
contrast, voltage pulses, which change the orientation, must be continually ap- 
plied. There is considerable interest in developing bistable cells, which support 
two (and possibly more) stable configurations with contrasting optical proper- 
ties. In bistable cells, power is needed only to switch between configurations. 
One mechanism for engendering bistability is to introduce microstructures into 
the geometry [71 [HI [22] . Nematic liquid crystals in cells with polyhedral features 
(eg, ridges, posts, wells) have been found to support multiple configurations. 
One such device, the PABN, or post-aligned bistable nematic cell, is shown in 
Fig. [T]([5]). It consists of a liquid crystal layer sandwiched between two planar 
substrates, with the lower substrate featured by an array of microscopic posts. 

As a simple model for such systems, we consider nematic liquid crystals 
in a convex polyhedron P C K"^ with orientation described by a director field, 
n : P — > RP^, taking values in the real projective plane. We consider the case of 
strong azimuthal anchoring, described by tangent boundary conditions. Tangent 
boundary conditions require that, on a face of P, n lies tangent to the face, but 
is otherwise unconstrained. It follows that on the edges of P, n is parallel to 
the edges, and therefore is necessarily discontinuous at the vertices. We are 
interested as to whether equilibria can be classified according to homotopy, and 
therefore restrict our attention to director fields which are continuous away 
from the vertices. For these, we can unambiguously assign an orientation to the 
director field (as P is simply connected), and regard n as a unit-vector field. We 
let Cj.(P, S*^) denote the space of continuous unit- vector fields on P satisfying 
tangent boundary conditions, or tangent unit-vector fields for short. 
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The elastic or Oseen-Frank energy of a configuration n is given by 



E 



JP 



[ [ifi(V ■ nf + K2{n ■ (V x n))^ + Ksin x (V x n))^ 



+ ■ {{n ■ V)n - {V ■ n)n)] dV. (1) 



Tangent boundary conditions imply that the contribution from the i4r4-term, 
which is a pure divergence, vanishes. We shall make use of the so-called one- 
constant approximation, in which the remaining elastic constants Ki, K2 and 
are taken to be the same and set to unity. In this case, ([T]) becomes the 
Dirichlet energy. 



Minimisers of the Dirichlet energy, which correspond to equilibrium configura- 
tions, are 5'^-valued harmonic maps, as are local minimisers, which correspond 
to metastable configurations. 

The homotopy classes of Cj^{P, S'^) are described in Section [2l and a lower 
bound for the infimum of the Dirichlet energy in each homotopy class is given 
in Section [31 The lower bound is expressed as a sum of minimal connections 
between fractional defects at the vertices of P, in analogy with the well-known 
result of [2] for the infimum Dirichlet energy of a set of point defects in M"^. For 
nonconformal homotopy classes, this bound can be improved by incorporating 
certain nonabelian homotopy invariants; this is shown explicitly for certain ho- 
motopy classes in a rectangular prism in Section 31 Unlike the case of point 
defects in R'^, the lower bound of Section O is expected to be strictly less than 
the infimum; achieving the lower bound would require concentration along a 
minimal connection, which would be incompatible with tangent boundary con- 
ditions. However, for P a rectangular prism, we can construct trial configura- 
tions in each homotopy class whose energies differ from the lower bound by a 
factor which is bounded independently of h (Section [5]). Generalising the con- 
struction to arbitrary P requires finding conformal maps on which preserve 
a given set of geodesies. 

It is an open question as to whether the infimum is achieved in a given 
homotopy class, as is the regularity of the local minimisers. Numerical re- 
sults presented in Section [51 suggest that some homotopy classes always contain 
smooth minimisers, while others may or may not depending on the geometry of 
P. Numerical results for a model of a bistable display suggests that the observed 
nematic configurations are topologically distinct. 

In addition to existence and regularity questions, it would be interesting 
to investigate namics under the influence of applied fields. Switching between 
configurations of different homotopy type requires the creation and destruction 
of defects, and one would like to understand this process in detail. 




(2) 
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Figure 2: Homotopy invariants, (a) Edge orientation (b) Kink number, n 
describes a 3/4-turn about the vertex (the image on is also shown), cor- 
responding to kink number —1. (c) Trapped area. The image of the cleaved 
surface on has signed area — 37r/2. 



2 Homotopy classification 

Given n e Cj^{P,S^) we can identify a number of discrete- valued quantities 
which depend continuously on n and which are therefore homotopy invariants. 
(Details may be found in [18] and [16j). Along an edge of P, n must lie parallel 
to the edge, so its value there is determined up to a sign, which we call an 
edge orientation (see Fig. [D^a)). Next, along a path on a face of P between 
two edges, n must lie tangent to the face, and therefore describes a geodesic on 
S"^, ie an arc of a great circle (see Fig. [21[b)). As the endpoints of the geodesic 
are fixed by the edge orientations, the geodesic may be assigned an integer- 
valued relative winding number, or kink number. By convention, the shortest 
geodesic is assigned kink number zero. Another invariant is associated with a 
surface which separates one of the vertices of P from the other vertices - we call 
this a cleaved surface (see Fig. [5l^c)). Along the boundary of a cleaved surface, 
n is determined up to homotopy by its edge orientations and kink numbers. 
Therefore, the signed area of the image of the cleaved surface itself, called the 
trapped area at the vertex, is determined up to an integer multiple of 47r (ie, 
some number of whole coverings of the sphere). 

Collectively, the edge orientations, kink numbers and trapped areas consti- 
tute a complete set of homotopy invariants for C^{P, S^); two configurations 
are homotopic if and only if their invariants are the same. We note that the 
invariants are not all independent - continuity of configurations on the faces of 
P implies that the kink numbers on each face satisfy a sum rule, while continu- 
ity on the interior of P implies that the trapped areas add up to zero. One can 
show that every set of invariants satisfying these sum rules can be realised. 

From the preceding discussion, it is evident that the invariants of n can be 
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Figure 3: Sectors for a tetrahedron. The great chcles of dhections tangent to 
the four faces partition the two-sphere into 14 sectors. 



determined from its values on a set of cleaved surfaces (the values of n on the 
corners and edges of the cleaved surfaces determine its edge orientations and 
kink numbers). Given a set of cleaved surfaces we can define an alternative set 
of invariants, the wrapping numbers, which will be used in subsequent sections. 
Let T C denote the set of directions which are tangent to one of the faces 
of P. Then T consists of a union of great circles. ~ T consists of a union 
of disjoint open spherical polygons, which we call sectors (see Fig. Let C° 
denote a cleaved surface separating the ath vertex of P, say, from the others, 
and let n° denote the restriction of n to C". Let E'^ denote the crth sector of 
S"^ . The wrapping number w""" is the number of times n° covers S"', counted 
with orientation. For n differentiable, this is given by 

= n^ix"^). (3) 

where lo is the area two- form on S"^ , normalised to have integral 47r, is the 
characteristic function of E'^, n* denotes the pull-back, and A'^ — |/g2 x'^ ^\ is 
the area of S'^. Alternatively, w""' can be expressed as the index of a regular 
value s G E"', ie 

w"" = sgn det(n'')'(r). (4) 

r I n'* (r)— s 

The wrapping numbers are homotopy invariants, and using Stokes' theorem 
can be expressed in terms of the edge orientations, kink numbers and trapped 
areas. These relations can also be inverted to obtain the edge orientations, 
kink numbers and trapped areas in terms of the wrapping numbers. Thus, the 
wrapping numbers constitute a complete (though redundant) set of homotopy 
invariants. 

If the nonzero wrapping numbers at a given vertex are all negative, the 
homotopy class is said to be conformal with respect to that vertex, and if pos- 
itive, anticonformal with respect to that vertex. A homotopy class is called 
nonconformal if there are vertices with wrapping numbers of different signs. 
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It is straightforward to count the number of invariants as well as the relations 
among them. Suppose that P has / faces, e edges and v vertices (so that, from 
Euler's formula, f—v+e ~ 2). Then P has v trapped areas, which satisfy a single 
sum rule; 2e kink numbers, which satisfy / sum rules; and e edge orientations. 
It also follows from Euler's formula applied to the set of tangent directions T, 
regarded as a graph on S"^, that there are generically (and at most) p — J + 2 
sectors ('generically' means that no direction is tangent to three or more faces of 
P). Thus, there are generically (and at most) {p — f + 2)v wrapping numbers, 
many more than the number of trapped areas and kink numbers. Among the 
constraints on the wrapping numbers, we point out that for a fixed sector E'^, 
their sum over vertices must vanish, ie 

a 

We will use h to denote both an admissible set of values of the invariants as well 
as the homotopy class in C^(P, S^) characterised by these values. 



3 Lower bound: minimal connection 

[5] established the infimum Dirichlet energy for unit-vector fields on M.^ with 
point defects of specified position and degree. The result is expressed in terms 
of a minimal connection between the defects, defined below. A similar argument 
yields a lower bound for the infimum Dirichlet energy for tangent unit-vector 
fields on P of fixed homotopy type, in which the vertices of P play the role of 
defects, and the wrapping numbers that of generalised degrees. Details may be 
found in [H], [B], and [IS]. 

We first review the result of [2j. Let = M'' — {r^, . . . , r"}, and let n : ^ 
denote a unit- vector field on fl. Continuous unit-vector fields on fl may be 
classified up to homotopy by their degrees, d = (d^, . . . , d") S Z", on spheres 
about each of the excluded points (the restriction of n to such a sphere may 
be regarded as a map from S'^ into itself). For n smooth, 

^— I n*uo (6) 

47r J|r-r3|=e 

for small enough e. For Vn square-integrable, the Dirichlet energy is given as 
in © by 

£;(n) = / {VnfdV. (7) 
Jn 

In order for -E'(n) to be finite, we require that 

^d^^O. (8) 

j 

Let C^(d) denote the homotopy class of continuous unit-vector fields with de- 
grees d satisfying ([5]), and let 

EfiHd) = inf E{n) (9) 

necO(d)nffi(o,S2) 
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denote the infimum energy in C^{d). 

Given two m-tuples of points in M^^ -p ^ (ai , . . . , a™) and TV = (b^ , . . . , b™) 
(whose points need not be distinct), a connection is a pairing (a^b"(j)) of points 
in V and Af, specified here in terms of a permutation tt G S,n {Sm denotes 
the symmetric group). The length of a connection is the sum of the distances 
between the paired points, and a minimal connection is a connection of minimum 
length. Let 

m 

L{V,Af) = min ^ ~ b'^^J) | (10) 

denote the length of a minimal connection, and let \d\ — ^ \d^. 

Theorem 3.1 [2] The infimum E^^^d) of the Dirichlet energy of continuous 
unit-vector fields on the domain fl = M.^ — {r^, . . . , r"} of degrees d^ about the 
excluded points is given by 

E^f{d)=^i^L{V{d),M{d)), (11) 

where V{d) is the \d\-tuple of excluded points of positive degree, with r'-' included 
d^ times, andN{d) is the \d\-tuple of excluded points of negative degree, with r'^ 
included \d''\ times. 

In fact, the result of [2] applies to more general domains with holes. 

Here we sketch the argument that 8TrL{V{d),J\f{d)) is a lower bound for 
E^f^^{d). It suffices to consider smooth unit- vector fields on fi, as these are dense 
in C^{d)r\H^{fl, S'^). For any orthonormal frame u, v, w, one has the inequality 

(Vn)2 > 2|(dCAn*tj)(u,v,w)| > 2(d^ A n*tj)(u, v, w), (12) 

where ^ is differentiable and \d£,\ < 1. p2|) follows from the fact that, at every 
point, there is at least one direction (say u) in which the directional derivative 
V„n := (u ■ V)n vanishes, while 

(V„n)2 + (V„,n)2 > 2|V,„n| \\/^n\ > 2|V„n A \/^n\. (13) 

Since = 0, it follows that c?^ A n*(jj ~ d{^n*LLj), so that 

(Vn)2 > 2d{^n*uj). (14) 

Substituting ((TJ]) into ([7]) and applying Stokes' theorem, we get a lower bound 

^;j?V)>2^erf^ (15) 

3 

which depends only on the values := ^(r-') of ^ at the defects. Since \d^\ < 1, 
these values are constrained by — £,'^1 < |r-' — r'^j. In fact, every set of ^■''s 
satisfying these constraints can be realised by a piecewise-differentiable function 
^ (eg, let ^(r) = maxj(^-' — |r — r-'D). Thus, one obtains a bound 



ET(d) > 2 max) £^d^, where > 0, - C\ < -r*^!, (16) 
z — ^ 

3 
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in the form of a finite-dimensional linear optimisation problem. 
The dual formulation is given by 

E'^^id) > 2min Vr^jfclr^' - r'=|, where r]jk > 0, V(r/jfe - ryfej) > ■ (17) 

jk k 

This is a sort of transport problem, in which the degrees are the quantities to 
be transported and the costs are the distances between defects. We can take 
rjjk to be unless > and d'' < 0. Without loss of generality, we may also 
assume that the degrees are either +1 or —1, so that there are an equal number, 
m := n/2, of each, with positions (a^, . . . , a™) and (b^, . . . , b™) respectively (if 
not, repeat each defect according to its multiplicity). (|17p becomes 

m 

E'^^id) > 2min ^ Mpq\aP - b«|, where Mpg > 0, ^ Mp, = ^ M^p = 1. 

p.,q^l p q 

(18) 

As M is constrained to be doubly stochastic, a theorem of Birkhoff pj implies 
that it lies in the convex hull of the set of m-dimensional permutation matrices. 
The optimal solution will be amongst the permutation matrices themselves, 
leading to E'^^{d) > BTTL{V{d),M[d)). 

A similar argument leads to a lower bound for the infimum Dirichlct energy 
for tangent unit-vector fields on P of given homotopy type. 

Theorem 3.2 [16j Let h = {w°'"} be an admissible topology for continuous tan- 
gent unit-vector fields on a polyhedron P. The infimum E^"-^{h) of the Dirichlct 
energy of continuous tangent unit-vector fields on P with invariants h is bounded 
below by 

Ep^{h) > 2A''L{r''{h),J\f^h)), (19) 

(7 

where (resp. ) contains the vertices of P for which is positive (resp. neg- 
ative), each such vertex included with multiplicity \w°''^\. 

Thus, to each sector a may be associated a constellation of point defects at the 
vertices v° with degrees The lower bound of (fT9|) is a sum of the lengths of 
minimal connections for these constellations, weighted by the sector areas A'^ . 

4 Lower bound: nonabelian invariants 

For nonconformal homotopy classes, the lower bound of Theorem 13.21 can be 
improved by incorporating certain nonabelian invariants. These invariants, and 
the sense in which they are nonabelian, are introduced in Section [4. II in a two- 
dimensional setting. For tangent unit-vector fields on P C M'' wc describe this 
phenomenon in a particular case (Section 14. 2p . reflection-symmetric homotopy 
classes in a rectangular prism. Details will be given in |15j . 
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4.1 Absolute degree and spelling length 

Let (j> • ^ be a smooth map of the two-disk into the plane. We recall 
that X € is a regular point oi (f> ii x is in the interior of and det (f>'{x) ^ 0, 
y e is regular value of (f> if all of its preimages are regular points, and a 
regular value has a finite number of preimages. Let TZ{(j)) C denote the set 
of regular values of (f>- From Sard's theorem, Imc/) — TZ{(j)) is of zero measure. 
Given y e 7^(</'), the algebraic degree of y (or degree, for short) is given by 

dAy)= sgndet(/.'(x), (20) 

and is invariant under smooth deformations of which preserve the boundary 
map d4>. We define the absolute degree of y by 

Clearly D^{y) is not invariant under all deformations which preserve dcj), and 

D^y) > my)\. (22) 

Let R = {yi, . . . , y„} C ^.{(j)) denote a set of n regular values of cj). We may 
regard the boundary dcj) as a map dcj) : M."^ — R from the circle to the 

n-times-punctured plane. Let 7ri(M^ — R,p) denote the fundamental group of 

— R, based at a point p, and let [d(j>] E 7ri(R^ — R,p) denote the homotopy 
class of d(f>. 

The fundamental group tti (M^ — R, p) may be identified with the free group 
on n generators, F{ci, . . . , c„) (see, eg, ilO]). Let us specify that the generator 
Cj corresponds to a loop which encircles yj once with positive orientation but 
contains no other y^'s. This determines the Cj's up to conjugacy. Given g g 
7ri(R^ — R,p), we define a spelling to be a factorisation of g into a product of 
conjugated generators and inverse generators, eg 

g ^ hlc'^h^^ ■ ■ ■ hr4';K\ (23) 

where hg G 7ri(R^ — R,p) and = ±1. The length of a spelling is the number 
of factors (ie, r in Define the spelling length, denoted An{g), to be the 

shortest possible length of a speUing of g (eg, A„(ciC2C^^C2"^) — 2). The spelhng 
length of the identity, e, is taken to be zero. It turns out that the spelling length 
of [dcj)] gives a lower bound on the sum of the absolute degrees of points in R. 

Proposition 4.1 Given cj) : R^ smooth, R = {j/i, . . . ,y„} C Ti{(j)), and 

7ri(R^ — R,p) ^ Fn{ci, . . . , On), with generators Cj as above. Then 



n 

Y.D^{y,)>Km]). 



(24) 



Let F„(ci, . . . , c„) denote the abehanisation of F„(ci, . . . , c„), obtained by 
taking all of the Cj's to commute, and given t; G F„(ci, . . . , c„), let ^ denote the 
corresponding element of Fn{ci, . . . , c„). Then g can be written as c^^ ■ ■ ■ c*" 
for some integers 6j. Let A„(g) = ■ Clearly A„(g) > An{g) (eg, 

A„(ciC2cJ"^c^^) = 0). It is readily seen that 

n 

J2\dM\=An{m). (25) 

Thus, Proposition l4.1l impHes that J^j ^tpiVj) is strictly greater than J2j \d<p{yj)\ 
provided that An{[d(j)]) is strictly greater than A„([90]). For example, if [d(j)] = 
ciC2cj~^c^^, then takes values yi or 2/2 at least twice, even though yi and 2/2 
are of degree zero. 

For our applications we shall want to consider maps v : — *■ 5^ from the 
two-disk into the two-sphere. Let R = {eo,...,e„} G Tlii^) denote a set of 
n + 1 regular values of v. In contrast to the case of maps to the plane, the 
algebraic degrees d{ej) are not determined by dv, since the image of v itself is 
determined only up to whole coverings of S^. We can remove this ambiguity 
by specifying the degree at one of the regular values, eg djy(eo) = do. We may 
identify the fundamental group 7ri(5'^ — R, q) with the free group F„(ci, . . . , c„) 
on n generators. As above, we specify that the generator Cj corresponds to a 
closed loop which encircles once with positive orientation but contains no 
other Bfc's, which determines the cj's up to conjugacy. cq, which corresponds to 
a loop about eo, may be expressed as a product of the generators ci through c„ 
and their inverses. In what follows, we write 6 ~ c to denote that b and c are 
conjugate. In analogy with Proposition 14. li we have the following: 

Proposition 4.2 Given u : — » 5^ smooth, R — {eo, . . . ,e„} C TZii^), and 

7ri(M^ — R,p) i2 Fn{ci, . . . ,Cn) with generators cj as above, such that cq G 
Fn{ci, . . . , c„). Suppose that (i^(eo) = rfo- Then 

n 

EDi^{ej) > \d()\ + min An{[du]gi ■ ■ ■ gr+\do\hi ■ ■ ■ K). (26) 

■'=° hi,...M,^r^Cg^ 

While it is straightforward to compute the spelling length of a given element g, 
evaluating (j26p may not be as straightforward. 

4.2 Reflection-symmetric homotopy classes in a prism 

A crude way to obtain a lower bound for the Dirichlet energy of tangent unit- 
vector fields on P is to estimate the contributions from nonoverlapping balls 
centred on each vertex. Let n denote a smooth tangent unit- vector field on P 
with invariants h. Let v° denote the ath vertex of P, and let O'^ C S'^ denote 
the set of directions about v" subtended by P. For r less than the length of any 
of the edges coincident at v°, define : 0° ^ by 

<(e) =n(v'^+re). (27) 
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Up to parameterisation, describes the restriction of n to a spherical cleaved 
surface of radius r. We have that 

|VnK+re)p > ^|K)'(e)|2 > 2^1 det«)'(e)|, (28) 
where the last inequality follows from the same reasoning as in ([T2|). Then 

E{n)>S2[ W^dr, where W° = / | det(z^^)' | de (29) 
a Jo Jo-^ 

and the RaS are chosen so that Ra + Rb < |va — v;,|. 

The quantity is just the unsigned area of Imi/°. The unsigned area of 
Im i/'^ n S"" is at least the area of 1]°' times the minimal absolute degree of the 
regular values in E'^. Thus we have that 

VF">V min D^a(e)A^. (30) 

^ ee7Z(i/»)nE'' 
a 

Noting that D^»(eo) > \w'"^° \ for aU Gq G TliK) ^ ^ we may apply Proposi- 
tion [i^] to (13DD to obtain 

> + y min A,^i{[du]gi ■ ■ ■ g,+|^a<,o|/ii ■ ■ ■ K) A" 

''=^"0 hx,...,hr~C-^ 

(31) 

(s in (pij) is the number of sectors). We note that it follows from ([50]) and (|22p 
that 

W^;" > y Iw;'"'!^''- (32) 

For certain homotopy classes of tangent unit-vector fields on a rectangular 
prism, i?, one can show that the estimate ([5^ based on spelling lengths leads 
to an improvement of the lower bound of Theorem 13.21 Let 

i? = {r |G < < Lj, j = a;,y,z}, (33) 

where for convenience we have chosen coordinates with the origin at one of the 
vertices and axes parallel to the edges. By convention, we take Lx > Ly > Lz- 
In this case, the sectors are the coordinate octants of with area A'^ = 7r/2. 

Reflection- symmetric homotopy classes on R are the homotopy classes of 
tangent unit-vector fields which are invariant under reflections through the mid- 
planes of the prism, 

n(a;, y, z) = n{Lx - x, y, z) = n{x, Ly - y, z) = n{x, y, L^- z). (34) 

In this case, the wrapping numbers at two vertices a and a related by a single 
reflection differ by a sign; 

w"" = -uf". (35) 
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Thus, the wrapping numbers about the origin determine all the rest, and for 
simplicity we denote these by w"^ . The prism, and reflection-symmetric config- 
urations in particular, will also feature in Sections O and [HI 

To estimate E^^\h) for reflection-symmetric ft, it suffices to consider tangent 
unit-vector fields n which are themselves reflection symmetric. From (j29p and 
([5^ it follows that 

>47r^|u;-|L., (36) 

CT 

which coincides with the lower bound (|19p of Theorem l3.2l (a minimal connection 
in this case is obtained by pairing vertices at the endpoints of the (shortest) Lz- 
edges of R). However, by using the estimate ((3T|) instead of ((32)) . we get the 
following. 

Theorem 4.3 |15j Let h he a nonconformal reflection- symmetric homotopy 
class in R. Let cr"*" denote the sector with largest positive wrapping number, 
denoted , and let a~ denote the sector with largest (in magnitude) negative 
wrapping number, denoted W~ . Let 

A(ft) = max ^ - x, \W-\^ ^ k"|-X, o|, 

y aeadj((7+) \w''>a ae adj {ij-} [w" <0 

(37) 

where adj(a) denotes the set of (three) octants adjacent to (ie, sharing an edge 
with) a, and x equal to or 1 depending on the signs of the edge orientations 
and kink numbers. Then 

Er^H) > A-K Iw"! + 2A(/i) j L,. (38) 
For typical nonconformal homotopy classes, A(ft) > 0. 



5 Upper bound in a prism 

In Theorem 13. 11 one obtains an equality for the infimum Dirichlet energy for a 
prescribed set of point defects, rather than just a lower bound, by constructing 
a sequence n^^-' whose energies approach 8T:L{V{d),J\f{d)). It can be shown that 
a subsequence n^*^^ approaches a constant away from lines joining the paired de- 
fects in a minimal connection (here assumed unique), while | Vn^*^-* p approaches 
a singular measure supported on these lines [2]. For tangent unit-vector fields 
on P, the boundary conditions preclude such a construction; n is required to 
vary across the faces of P, and therefore throughout its interior. However, by 
constructing tangent unit-vector fields which saturate the local inequality (|12p 
over most of P, we can produce upper bounds for the Dirichlet energy with 
the same scaling with homotopy invariants as the lower bound of Theorem 13.21 
Details are given in [T31 El [ini [H] . 
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Here in outline is a procedure for constructing such configurations. Fix a set 
of values h of the homotopy invariants. As in Section [42l let C 5^ denote 
the set of directions about the vertex v° subtended by P. Define spherical 
cleaved surfaces 

C = {w" +r''e\ee S^}, (39) 

where r" is taken to be less than half the length of the smallest edge coincident 
at v° (so that the C""s do not intersect). Specify n on C° so as to satisfy 
tangent boundary conditions with wrapping numbers given by h, and take n to 
be constant along rays from C" to v". It remains to define n on P, the (closed) 
domain obtained by excising the cones between the C°'s and v^'s. The boundary 
of P is composed of i) the C°'s and ii) the faces of P truncated by the C""s. 
Extend n smoothly to these truncated faces so as to satisfy tangent boundary 
conditions. Choose a point p in the interior of P. Along rays from C° to p, take 
n to be constant. Along rays from each truncated face to p, rotate the values of 
n out of the tangent plane to the outward normal. There emerges a discontinuity 
at p, but this is easily removed. If n is specified on the C""s to be conformal or 
anticonformal, except possibly on a small subset where its derivative is suitably 
controlled, then the local inequality is saturated throughout most of P, 
and the Dirichlet energy can be shown to be proportionate to the lower bound 
of Theorem 1321 independently of h. 

The main difficulty in carrying out this procedure is in defining n on the C""s. 
Let u"^ : O"^ ^ S'^ be given by /^"(e) = n(v'' + r'^e) (similarly to ^7^). We 
note that O" is a geodesic polygon on S'^; its sides are arcs of the great circles 
of directions tangent to the faces of P which are coincident at v". Tangent 
boundary conditions require that v'^ maps each side of 0° into the great circle 
containing it. If /i is conformal with respect to (the anticonformal and 
nonconformal cases are discussed below), we are led to the following: 

Problem 5.1 Find conformal maps on which preserve a given set of geodesies. 

Restricting the domain of such a map to O'^ yields a candidate for i/". 

In the case of the rectangular prism R, Problcm l5.1l is readily solved. There 
are three geodesies which meet at right angles, and which may be taken to 
be the great circles about x, y, and z. Under the stcreographic projection 
e I— > w = (cx + iCy) / {1 + Cz) from to the extended complex plane, these are 
mapped to the real axis, imaginary axis and unit circle respectively. Problem l5.ll 
becomes one of finding locally analytic functions f{w) such that i) f(w) is real 
when w is real, ii) f{w) is imaginary when w is imaginary, and iii) \f{w)\ = 1 
when \w\ = 1. Property i) implies that / is real; ii) then implies that / is odd; 
iii) then implies that f{l/w) = 1/ f{w). Therefore, if is a zero of /, then —w^ 
and iD* are zeros, while l/w* is a pole. Restricting to / to be meromorphic, we 
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may conclude that / is rational of the form 



fH 




n( (-^-^')(-^-^^) V (40) 



X 

1=1 



The Tj's denote the real zeros {pj = 1) and poles {pj — —1) of / between and 
1; the Sfe's, the imaginary zeros and poles of / (according to cTfc = ±1) between 
and i ; and the ti^s, the complex zeros and poles of / (according to r; = ±1) 
with modulus less than one and argument between and tt/2. 

The parameters in ([^0]) can be chosen to realise any admissible set of con- 
formal (ie, nonpositive) wrapping numbers. Anticonformal topologies can be 
realised by replacing w with w. Nonconformal topologies can be produced by 
modifying / in a small neighbourhood to be anticonformal and smoothly inter- 
polating between the conformal and anticonformal domains. 

Let 

Epih) - ^2A-L(p-(/i),A/-'^(M) (41) 

(T 

denote the lower bound of Theorem 13.21 

Theorem 5.2 p,6j Let R denote a rectangular prism with sides of length > 
Ly > Lz and largest aspect ratio k = L^/Lz- Then 

E'^f{h) <CK^Ep{h) (42) 

for some constant C independent of h and L^, Ly, Lz- 

In the proof of Theorem l5.2i the positions of the zeros and poles of the conformal 
map (|40p must be chosen carefully to ensure that the bound is achieved. 

For reflection-symmetric conformal homotopy classes (cf (|34p). a simpler 
construction leads to an improved result, in which C = 1 and is replaced by 
{Ll + Ll + Llf/^/Lz. 

Theorem 5.3 T3] Let R denote a rectangular prism with sides of length L^ > 
Ly > Lz and h a reflection- symmetric homotopy class which is conformal about 
one of the vertices. Then 



(t2 , t2 , 7-2^1/2 

gjT^W < ^ y Epih). (43) 

L, 



Theorem 15.31 extends to nonconformal reflection-symmetric homotopy classes, 
provided Ep{h) is replaced by the lower bound given by Theorem [ 
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Figure 4: Unwrapped configuration in class hg 

6 Existence and regularity of local minimisers: 
numerical results 

Using direct methods, one might expect to estabhsh the existence of a global 
minimiser of the Dirichlet energy for tangent unit-vector fields on P. The exis- 
tence of (continuous) local minimisers in a given homotopy class h, however, is 
more difficult to address. The homotopy classes are not weakly closed, so that 
the existence of such local minimisers is not guaranteed; Ep^h) may not be 
realised (just as the infimum energy for a prescribed set of point defects in R'^ is 
not realised). We also recall the Hardt-Lin phenomenon ^6] - global minimisers 
of the Dirichlet energy may have interior singularities, even when continuous 
unit-vector fields are admissible. If continuous local minimisers exist, then one 
would like to analyse their regularity [191 UHl HI HI] • 

Questions about the existence and regularity of continuous local minimisers 
of given homotopy type appear to be open for the problems we are considering. 
Below we describe some numerical results which suggest that, for some homo- 
topy classes, smooth minimisers always exist, while for others, they may exist 
or not depending on the geometry of P. 

The first examples concern two reflection-symmetric homotopy classes in a 
rectangular prism, denoted here by Hq and hi, which are both conformal with 
respect to one of the vertices (and, therefore, conformal or anticonformal with 
respect to the others). Details are given in 13J. ho is the simplest possible, 
in which there is a single nonzero wrapping number equal to —1, so that n 
takes values in a single octant of S^. The restriction of n to a spherical cleaved 
surface corresponds to the conformal map given by /o(w) = w (cf (jM])). Such 
a configuration is shown in Fig. [H The lower bound for the infimum energy of 
Theorem 13.21 is -inLz- The upper bound of Theorem 15.31 can be improved in 
this case by explicit evaluation of the Dirichlet energy for a trial configuration, 
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yielding 




1 1 1 1 1 y 

' 2 ' 2 ' 2 ' 2 ' 7-2 




) 



+ {x 



z) + {x ^ z ^ y), 



(44) 



where F2{a, (3, /?', 7, 7'; s, i) is the Appell hypergeometric function [S]. For a unit 
cube, we get the bounds 



We computed minimisers numerically using two methods, namely solution of 
the Euler-Lagrange equation (using FEMLAB, a commercial PDF solver) and 
gradient descent. The converged energies from both methods agree, giving ap- 
proximately 14.8. The converged unit- vector field is indistinguishable from Fig. 
ID at the resolution shown, and appears to be regular away from the vertices. 

The homotopy class hi is the next simplest among the reflection-symmetric 
conformal classes. There are three nonzero wrapping numbers equal to —1 in 
contiguous octants, so that n takes values in three-quarters of a hemisphere. 
The restriction of n to a spherical cleaved surface corresponds to the conformal 
map 



Such a configuration is shown in Fig. [5] On the xy-faces of the prism, n executes 
a three-quarter turn about each vertex, corresponding to kink numbers of ±1; 
the kink numbers on the other faces all vanish. As the parameter s approaches 
1 from below, half-turns becomes concentrated along the y-edges, while /i ap- 
proaches /o away from the y-edges. The corresponding family of configurations 
n^. is weakly but not strongly continuous with respect to s (an example of the 
fact that hi is not weakly closed). 

Both numerical methods indicate that hi supports a smooth local minimiser 
for sufficiently thin slabs [Lx/Ly^Lx/Lz < 1/10), while for aspect ratios closer 
to unity, the numerical solution converges to the minimiser in /iq. Some insight 
into this behaviour is provided by computing the Dirichlet energy of trial con- 
figurations characterised by the one-parameter family (j46|) , as shown in Fig. [6l 
For a cube (dashed curve), the energy approaches a minimum as s approaches 
1, corresponding to a configuration in which half-turns concentrate along the 
y-edges. For Lx — 20, Ly — 10, = 1 (solid curve), the energy has a minimum 
for s between and 1, corresponding to a smooth configuration. Note that 
concentration along the shortest [LzYedges (which support the minimal con- 
nection) is not compatible with the topology, as the nonzero kink numbers lie in 
the xy-faces. Analogous arguments suggest that refiection-symmetric homotopy 
classes with two or more nonzero kink numbers do not contain smooth minimis- 
ers (for these classes, concentration along the shortest edge is compatible with 
the topology). However, it is conceivable that more non-reflection-symmetric 



12.5 < E'^^ha) < 15.3 



(45) 
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Figure 5: A configuration n in hi, generated by the conformal map /(w) = 
w{w'^ + s^)/(s^«;^ + 1) for s = .5. n describes a 3/4-turn about each vertex in 
the top and bottom face. 




Figure 6: Scaled energy e = E/{LxLyLzY^^ of the conformal configuration 
w{w'^ + s'^)/{s'^w'^ + 1). Solid curve: = 20, Ly = 10, = 1. Dashed curve: 
Lx = Ly = Lz = 1. 
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homotopy classes (for which minimal connections do not necessarily pair vertices 
along edges) support smooth local minimisers. 

The last numerical example is an idealised model of the PABN device. De- 
tails are given in [T^]. In fact, the model lies outside the class of problems we 
have considered so far; the domain is not a polyhedron, and the boundary con- 
ditions are not purely tangent. We take the PABN to consist of a rectangular 
post of square cross-section centred on the bottom surface of a rectangular cell 
of square cross-section, as in Fig. [T] In keeping with the device dimensions, the 
cell height is taken to be three times the cell width, and the cell width to be twice 
the post width. The height of the post is variable. Boundary conditions are 
dictated by material characteristics of the substrates. Tangent boundary con- 
ditions apply on the bottom substrate and on the post, while normal boundary 
conditions are appropriate for the top substrate. Periodic boundary conditions 
are imposed on the vertical sides of the cell, simulating a two-dimensional ar- 
ray of cells supporting the same nematic configuration (at a given time) and 
comprising a single pixel. 

We consider four simple homotopy classes, in which the kink numbers are 
zero and the trapped areas taken to have their minimal allowed values. The 
orientation of n on the horizontal edges of the post are fixed, as in Fig[7Ka). 
The classes are distinguished by the relative orientations of n on the vertical 
edges of the post. Up to symmetry, there are four distinct possibilities. For the 
tilted class T, the orientation on all four vertical edges is the same. The other 
three classes, called planar, are obtained by taking the orientation to be opposite 
on, respectively, one of the vertical edges (the Pi class), two adjacent vertical 
edges (P2), and two opposing vertical edges (Pa). Configurations in T exihibit 
a large vertical component Uz in the region around the post. In configurations 
in Pi - P3, riz is suppressed by the change in orientation between the vertical 
edges. 

Local minimisers for each of these homotopy classes were computed using 
FEMLAB for a range of post heights. The converged configurations appear to 
be smooth away from the vertices of the post. In Figure [5] we plot the converged 
energies of the local minimisers as a function of post height. The tilted class 
has the lowest energy, which is consistent with experimental observations which 
show that the liquid crystal always relaxes into the high-tilt state when cooled 
from the isotropic state ^ . The computations support the hypothesis that the 
bistable states of the PABN are topologically distinct. 
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Figure 7: Edge orientations for the four PABN configurations, (a) Orientations 
on the horizontal edges are the same for alL (b) Tilted profile T. n points up 
on all vertical edges, (c) Planar profile Pi. n points down on a single vertical 
edge (d) Planar profile P2- n points down on a pair of adjacent vertical edges 
(e) Planar profile P3. n points down on a pair of opposite vertical edges. 
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